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ABSTRACT 

The first hydrostatic core (FHSC) represents a very early phase in the low-mass star formation 
process, after collapse of the parent core has begun but before a true protostar has formed. This large 
(few AU), cool (100 K), pressure supported core of molecular hydrogen is expected from theory, but 
has yet to be observationally verified. Here we present observations of an excellent candidate for the 
FHSC phase: Per-Bolo 58, a dense core in Perseus that was previously believed to be starless. The 
70 /jm flux of 65 mJy, from new deep Spitzer MIPS observations, is consistent with that expected for 
the FHSC. A low signal-to-noise detection at 24 /im leaves open the possibility that Per-Bolo 58 could 
be a very low luminosity protostar, however. We utilize radiative transfer models to determine the 
best-fitting FHSC and protostar models to the spectral energy distribution and 2.9 mm visibilities of 
Per-Bolo 58. The source is consistent with a FHSC with some source of lower opacity through the 
envelope allowing 24 /im emission to escape; a small outflow cavity or a cavity in the envelope are 
both possible. While we are unable to rule out the presence of a protostar, if present it would be one 
of the lowest luminosity protostellar objects yet observed, with an internal luminosity of ~ O.OILq. 
Subject headings: stars: formation — infrared: ISM — submillimeter: ISM — radiative transfer 



1. INTRODUCTION 

One stage of the low mass star formation paradigm has 
yet to be observed: the first core or first hydrostatic core 
(FHSC). During this brief phase between the initiation 
of core collapse and the appearance of a "true" protostar, 
the growing central core of molecular hydrogen (H2) is 
heated to the point at which it obtains hydrostatic bal- 
ance. When the central temperature of the first core 
reaches approximately 2000 K, the H2 is dissociated and 
a second collapse forms the true protostar. 

Th e FHSC has long been expected from theory (|Larsonl 
1969T). but has n ot yet been observationally verified. 



Chen et "ail (|2010f ) recently presented a promising can- 



didate, L1448 IRS2E, based on detection of a CO out- 
flow without a corresponding mid-infrared source. The 
observed outflo w is faster than wou ld be expected for a 
typical FHSC (|Machida et a l. 2008), but might be con- 
sistent with the very end of the FHSC phase. 

The expected effective temperature of the FHSC, 
which has a radius of several AU, is approximately 100 K; 
as the majority of this cool emission is absorbed and re- 
radiated by the even cooler envelope, the emergent spec- 
trum of the FHSC is very similar to that of a starless 
core. The strongest observational signature is likely to 
be an increase in luminosit y for A ~ 40 — 100 /im as 
compared to starless cores (IBoss fc Yorkd 11991 . Fur- 
thermore, given the expected short lifetime of the FHSC 
(t - 10 3 -3 x 10 4 vr: IOmukail[2lffl7l IBoss fc Yorkdll99"l . 
there should only be a few observable even in large sam- 
ples of cores and protostars. Based on the relative life- 
times of the FHS C and true protostars (5.4 x 10 5 yr; 
lEvans et al.ll2009l ). we expect only one FHSC for every 
18 to 540 protostars. Before large Spitzer surveys such 



as the "From Molecular Cores to Planet-for ming Disks" 
Lega cy Program ("Cores to Disks" or c2d; lEvans et al.l 
2003) , the number of known protostars in any given re- 
gion was typically too low to expect even one FHSC. 

We have recently completed a census of starless 
cores and embedded protostars in the Perseus, Ser- 
pens, and Ophiuchus molecular clouds (jEnoch et al.1 
120081 l20"09af ) . bas ed on large-scale B o locam 1.1 mm con- 
tinuum surveys (Enoch et a l. 2006; lYoung et al.l 120061 : 
lEnoch et al.ll2007l) and IRAC and MIPS maps from c2d 
(iJprgensen et al.l 120061: iHarvev et alj|2006t iRebull et al.l 
l2007HHarvev et al.ll2007bHPadgett et al.H2008t ). Perseus. 
Serpens, and Ophiuchus each harbor between 30 and 70 
embedded protostars - objects after the second collapse, 
and thus more evolved than FHSCs, but still embedded 
in and accreting from their natal cores - and should con- 
tain at least a few FHSCs (if the lifetime is longer than 
10 4 yr). Spitzer was sufficiently sensitive to detect 70 /tin 
emission from a 0.015 M Q FHSC at a distance of 250 pc 
(Omukai 2007), but large surveys such as c2d were not 
deep enough at 70 /im to detect the FHSC. 

In this Letter, we present new deep Spitzer MIPS ob- 
servations that reveal a candidate FHSC in the Perseus 
molecular cloud. Utilizing literature data, we compare 
the spectral energy distribution (SED) and millimeter 
visibilities of this source to radiative transfer models of 
a FHSC and very low luminosity protostar. 

2. TARGET SELECTION 

Figure Q] shows the expected observed SEP of a FHSC, 
from the radiative transfer models of lYoung fc Evansl 
(2005). Based on these models, an unambiguous detec- 
tion of the FHSC requires detection of compact 70 /im 
emission in a dense core with no emission at A < 24 /im. 
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The 70 /im sensitivity limit of the c2d survey was ap- 
proximately 75 mjy (5a), somew hat higher than the 
flux expected for a ty pical FHSC (|Boss fc Yorkd 119951: 
I Young fc Eva^\2()m . Thus any FHSCs appeared as 
starless cores in our census, and deeper 70 urn follow-up 
observations were necessary to reveal themQ 

We selected targ ets from our census of starless cores 
(|Enoch et al.ll2008f ) . which includes 108 cores and is com- 
plete to M core > 0.2M©, based on three criteria. (1) Vi- 
sual inspection of the c2d 70 /im maps at starless core 
positions, looking for faint emission below the formal 
5(7 detection limit. (2) Starless cores with mass-to-size 
ratios above the critical value for stable Bonnor-Ebert 
spheres (central to surface density ratio p c /po > 14, or 
M [Mq]/FWHM [pc] > 31.9 for T = 10 K). These cores 
are likely to be collapsing, but the lack of near- or mid- 
infrared emission indicates that they have not yet formed 
a true protostar. (3) Spatially compact starless cores, ei- 
ther unresolved by the 31" Bolocam beam or compact in 
follow-up CARMA 2.9 mm interferometric maps (Enoch 
et al., in prep). Again, these cores have high central 
density and may be collapsing. 

In the three clouds, 7 starless cores met one or more 
of the above criteria. This is by no means a complete or 
unbiased sample; rather, we used these criteria to identify 
interesting objects for deeper follow-up observations. 

3. OBSERVATIONS AND RESULTS 

Deep 70 /im maps of our 7 targets were obtained with 
the Multi-Band Imaging Photo meter for Spitzer (M IPS) 
on the Spitzer Space Telescope (jWerner et al.H200l) dur- 
ing 2008 October 20 - October 27. Integration times 
ranged between 1100 — 2000 s pixel -1 , depending on the 
brightness of the background estimated from c2d maps. 
The resulting 5a detect ion limit of ~ 7 .5 mjy was chosen 
based on the models of lOmukail (|2007D . to detect FHSCs 
with mass > 0.015 M Q at the distance of Perseus and 
Serpens (250 - 260 pc). 

Da ta were reduced using the c2d pipeline (|Evans et al.l 
2007). One source was clearly detected at 70 /im: 
[EYG2006] Bolo 58 (RA=03 29 25.7, Dec=+31 28 16.3), 
hereafter Per-Bolo 58 , a den se co re identified as s tarles s 
by both lEnoch et al] (|2006[ ) and IHatchell et al] (|2007|) . 
Per-Bolo 58 was included in our target list based on faint 
emission in the c2d 70 /im map and a spatially com- 
pact 1.1 mm core. It is located in the outskirts of the 
NGC 1333 young cluster (see Figure [2 right). 

The 70 /im flux of Per-Bolo 58 is 65 ± 6 mjy. A subse- 
quent search of the c2d deep catalogs (jEvans et al.H2007|) 
revealed a low signal-to-noise 24 /im detection at the po- 
sition of Per-Bolo 58, with a flux of 0.88 ± 0.24 mjy. 
This source was not included in the "high reliability" 
or "YSO candidate" c2d catalogs, as it did not meet 
the 7a signal-to-noise requirement. Given the good po- 
sitional correspondence and point-like emission, we be- 
lieve the 24 /im source is associated with Per-Bolo 58. 
There is also a nearby low signal-to-noise 4.5 /im detec- 
tion (0.024 ± 0.007 mjy) that we do not believe is as- 
sociated, as the only point-like emission is offset by 10" 
from the 24 /im and 70 /im sources. We treat the 4.5 /im 
point as an upper limit. We estimate a 160 /im flux of 

1 We did an initial check of the c2d catalogs for starless cores 
detected only at 70 /im but found no such objects. 



2.8 Jy from the c2d MIPS 160 /im map; this value has a 
large uncertainty (~ 50%) due to non-uniform extended 
emission in NGC 1333. Spitzer images at 4.5, 24, 70, and 
160 /im are shown in Figure [2J 

All known continuum fluxes for Per-Bolo 58 are given 
in Table [TJ and the observed SED is shown in Fig- 
ure [T] In addition to the Spitzer data described above, 
we include SHARC II 350 /im (Dunham et al., in 
prep), SCUB A 850 urn. (IHatchell et al.l l200l. and Bolo- 
cam 1.1 mm (jEnoch et al.l I2006[) fluxes, a s well as the 
CARMA 2.9 mm map from ISchnee et al.l (pOlOlh The 
interferometric CARMA 2.9 mm map provides both a 
long wavelength flux measurement and a measure of the 
resolved radial intensity profile. The 2.9 mm visibility 
amplitudes versus uw-distance are shown in Figures [3] 
andgj 

4. DISCUSSION 

Our 70 /im detection of Per-Bolo 58 indicates the pres- 
ence of an internal luminosity source in this dense core, 
with a flux consistent with that expected for a FHSC. 
The detection at 24 /im is inconsistent with a FHSC sur- 
rounded by a spherically symmetric envelope, however 
(Figure [I| , leaving open the possibility that Per-Bolo 58 
may be a very low luminosity protostar. We discuss these 
possibilities below. 

4.1. Radiative Transfer Models 

We use RADMC, a two-dimensional Monte Carlo ra- 
diative transfer code (Dul lemond fc Dominikl l2004j) to 
model a FHSC and low luminosity protostar. We adopt 
a density profile that includes a rotating infalling enve- 
lope characterized by total mass M env , centrifugal r adius 
Rcent, and outer radius R out (e.g. lUlrich et al.fl967D . and 
outflow cavity characterized by full opening angle a. The 
envelope mass is held fixed at M env = 1.2 Mm, as cal- 
culated from the 1.1 mm flux (see lEnoch et all 1200(1 
assuming a dust temperature of T = 8 K, consistent 
with radiative transfer models of such a low luminosity 
sour ce. Details of the de nsity profile and modeling are 
as in lEnoch et all (|2009bf ). 

Dust opacities are from Table 1, column 5 of 
IQssenk opf & Henninal (|1994D for dust grai ns with thin ice 
mantl es, including scattering. Following IDunham et all 
(2010), we remove the effects of excess back-scattering 
of the interstellar radiation field (ISRF) by subtract- 
ing the flux for a star l ess co re model at short wave- 
lengths. IDunham et al.l ()2010D make this correction for 
A < 10 /im; here we additionally include the small 
plateau of ISRF emission at A = 10 — 40 /im. While 
this correction affects the shape of our model SEDs at 
short wavelengths, the correction at 24 /im is small, so it 
does not significantly affect the model fits. 

To match the peak of the SED at 80 /im < A < 
1000 /im we include an ISRF th at is a factor of 3 stronger 
than the "Black-D raine" ISRF (jEvans et al.ll200il: IBlackl 
119941 : IDraind ll"978). A more intense radiation field may 
be expected in NGC 1333, where there is significant ex- 
tended emission at A < 24 /im (e.g. Figure [5]). Increas- 
ing the ISRF improves the overall fit of all models, but 
does not change the best-fitting envelope parameters, be- 
cause the SED peak depends primarily on fixed parame- 
ters (M env , ISRF). 
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4.2. Per-Bolo 58 as a First Hydrostatic Core 

Our FHSC models include a cool (T eff = 100 K), large 
(R ~ 2 AU, where the actual radius is determined from 
Lint — ^R 2 o~T e f f A ) internal luminosity source. We run 
a grid of models varying L mt (0.006, 0.012, 0.025 L Q ), 
Rout (4000, 6000, 8000, 10000 AU), R cent (10, 50, 100, 
200, 300 AU), a (0, 5, 10, 20, 40, 60 deg), and inclination 
(i; 5, 10, 15. ..90 deg). 

Models are compared to the observed SED of Per- 
Bolo 58 using a x 2 analysis to constrain Li nt , R ce nt, 
a, and i. An internal luminosity of 0.012 L0 is clearly 
preferred by the SED, and we derive other parameters 
with Li n t fixed. The SED is rather insensitive to R ou t, 
so we also compute the visibility amplitudes of each 
model as a function of uv -distance, and compare to th e 
2.9 mm CARMA observations from lSchnee et all (|2010D . 
At small ww-distances (< lOfcA), the visibilit ies are sensi- 
tive to t he spatial extent of the envelope (see lEnoch et al.1 
(|2009bD for more details), and we find a best fit for 
R out = 8000 AU. 

The best fitting FHSC model (L mt = 0.012 L©, T eff = 
100 K, R = 1.7 AU, R cen t = 50 AU, R out = 8000 AU, 
a = deg, i = 25 deg) is compared to the observed SED 
and 2.9 mm visibilities in Figure [3l Another well-fitting 
model with R cen t = 10 AU and a = 40 deg is also shown. 
A number of models have similar x 2 values, but all have 
either an envelope cavity (Rcent > 50 AU) or an outflow 
(a > 20 deg). Models without an envelope cavity are 
clearly unable to match the observed 24 /im point, due 
to the high opacity through the envelope; such a model 
with R C ent = 10 AU and a = deg is shown in Figure |3] 
for comparison. 

Although an envelope cavity is not expected from sim- 
ple physical models of the FHSC, in which the envelope 
extends down to the FHSC and no outflow is present, 
there are several cases in which we might expect a cavity. 
MHD calculations have shown that the FHSC can drive a 
molecular outflow (e.g. iMachida et al"1l2008l ). Such early 
outflows should be weak, but could evacuate enough en- 
velope material to allow a small amount of 24 /im flux 
to escape. Similarly, a larger centrifugal radius might 
be expected if the inner envelope were cleared by a bi- 
nary source or weak outflow. We note that 4.5 /xm is the 
Spitzer band associated with shock emission in outflows, 
and the nearby weak 4.5 /jm emission could be tracing a 
weak outflow. 

4.3. Per-Bolo 58 as a Very Low Luminosity Protostar 

A number of dense cores previously believed to be star- 
less were found by Spitze r to harbor yery lo w luminos- 
ity objects (VeLLOs; e.g. lYoung et all 120041 ). The best 
current explanation of VeLLOs is that they are true pro- 
tostars, but with low internal luminosities due either to 
very low masses or very low accretion rates. Most VeL- 
LOs have internal lum inosities between 0.02 and 0.1 L Q 
(jDunham et al.ll2008lh although the 24 and 70 fxm fluxes 
of Per-Bolo 58 are lower than typical VeLLOs, it may be 
an extreme example of a low luminosity protostellar pop- 
ulation. At least two other protostellar sources with simi- 
lar luminosities are known: a candidate embedded proto- 
brown d warf in Taurus with L int ~ 0.003 L Q (SSTB213 
J041757; iBarrado et al.|[2009h . and Cha-MMS 1, which 
has an estimated internal luminosity of 0.01 — 0.02 L 



(|Belloche et alJl200"flPl 

To model Per-Bolo 58 as a protostar, we use a sim- 
ilar grid as above, but include a hot internal luminos- 
ity source typical of a true protostar (T e ff = 3000 K, 
R ~ 1 R©). The best fitting protostar model (Li n t = 
0.012 L , Rcent = 50 AU, R out = 8000 AU, a = 10 deg, 
i = 30 deg) is shown in Figure |U Protostar models 
with relatively high inclinations (i > a) are nearly indis- 
tinguishable from the corresponding FHSC models; in 
both cases essentially all of the emission from the central 
source is reprocessed by the envelope, masking the de- 
tails of the internal luminosity source. We can rule out 
protostar models with low inclinations, however, which 
produce far too much emission in the near-infrared. 

Again, an internal luminosity of 0.012 L Q is strongly 
favored. This value is corroborated by the relation- 
ship between inter nal luminosity and ob served 70 /an 
flux determined bv iDunham et all ([2008): L; n t = 3.3 x 
10 8 F 7 ° 94 Lq, where F 70 is the flux ( m erg cm s J 
at 140 pc. For Per-Bolo 58 this relationship yields 
Li n t = 0.014 Lq. If Per-Bolo 58 is a true protostar then 
it is one of the lowest luminosity embedded protostars 
known, with Li n t ~ 0.01 Lq. 

4.4. Possible Observational Tests 

We are currently unable to conclusively determine the 
evolutionary state of Per-Bolo 58. W hile the Herschel 
Gould Belt survey (An dre et alll2010l ) will further refine 
the SED, additional observations are needed to distin- 
guish a FHSC with a small envelope cavity from an ex- 
tremely low luminosity protostar. Possible observational 
tests include measuring the outflow velocity and finding 
evidence of water ice evaporation. 

Numerical MHD models suggest that the FHSC can 
drive a molecular o utflow, but with quite low veloc- 
ity (v ~ 3 km s" 1 ; IMachida et "all 120081 : iTomida et all 
120101 ). As protostell ar outflows typically h ave speeds of 
10 - 20 km s" 1 (e.g. lArce fc Sargentll2006l ). detecting an 
outflow and measuring its velocity, either in the millime- 
ter or with Herschel spectroscopy, could test our hypoth- 
esis that Per-Bolo 58 is a FHSC. No outflow was detected 
in the HARP 12 CO (3-2) study of lHatchell fc Dunham! 
(2009), but a FHSC outflow should be below their detec- 
tion criteria (linewing strength above 1.5 K at 3 km s _1 
from the core velocity). 

Unlike the FHSC, for which the maximum dust tem- 
perature in the envelope is only ~ 80 K, temperatures in 
our protostar models are high enough in the inner few AU 
(up t o 300 K) to evaporate water ice (e.g. iFraser et all 
2001). ALMA is capable of detecting the resulting H 2 
or HDO lines, around 300 GHz and 240 GHz, even with 
the beam dilution of an 0.08" (20 AU at 250 pc) beam. 
Given the outflow non-detection at the level expected for 
a typical protostar, finding chemical signatures indicative 
of high gas temperatures may be the best way to rule out 
Per-Bolo 58 as a FHSC. 

5. CONCLUSIONS 

In an effort to observationally verify the theoretically 
predicted FHSC phase, we have obtained deep Spitzer 

2 IBelloche etTaTI !|2009 ) note that Cha-MMS 1 may also be a 
FHSC candidate, however. 
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70 /Ltm images of a small sample of dense starless cores. 
Per-Bolo 58 is detected at the level expected for the 
FHSC, making it a very promising FHSC candidate. A 
weak detection at 24 /xm leaves open the possibility that 
it might be an extremely low luminosity protostar, how- 
ever. 

We are able to reproduce the observed SED and 2.9 
mm visibilities of Per-Bolo 58 with radiative transfer 
models of a FHSC, although some source of lower en- 
velope opacity allowing 24 ^m emission to escape is re- 
quired: cither a small outflow cavity or a spherical cavity 
in the envelope of ~ 50 AU. While we cannot rule out 
the possibility that Per-Bolo 58 has already formed a true 
protostar, with an internal luminosity Li nt ~ 0.01 L Q it 
would be one of the lowest luminosity protostars known. 
Additional observational tests to clarify the evolution- 
ary state of Per-Bolo 58 include measuring the velocity 
of any outflow present, and looking for evidence of wa- 
ter ice evaporation in the inner few AU. Finally, we note 
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TABLE 1 
Observed SED of Per-Bolo 58 



Wavelength 


Flux 


cr(Flux) 


Aperture 


Notes 


(H 


(mjy) 


(mjy) 


diameter (") 




4.5 


0.024 


0.07 


2.2 


upper limit; Spitzer (c2d) 


24 


0.88 


0.24 


7 


Spitzer (c2d) 


70 


65 


6 


17 


Spitzer; this work 


160 


2870 


1600 


40 


Spitzer (c2d) 


350 


6100 


1200 


10 


SHARC II; Dunham et al., in prep 


850 


920 


200 


18 


SCUBA; Hatchell et al. (2005) 


1100 


330 


30 


40 


Bolocam: Enoch etld~72006) 


2930 


13 


6 


15 


CARMA; Schnee et al. (2010) 



Note. — Continuum fluxes are calculated by aperture photometry within the aperture listed, with the ex ception of the 2.9 m m flux, 
which is determined by a Gaussian fit (FWHM= 15"). Our 2.9 mm flux does not include the SZA data from Schnee et al. (2010), which 
probes larger spatial scales than we include in our models. 
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Fig. 1. — Radiative transfer models from Young & Evans (2005), for a 1 Mq core at the distance of Perseus (d = 250 pc). SEDs 
represent a time sequence, from a starless core (dotted line) to FHSC (dashed lines) to true protostar (solid lines), with short wavelength 
flux increasing monotonically with time and time steps of 2000 yr. The observed SED of Per-Bolo 58 (diamonds; see Section O includes 
fluxes from Table Q] The 70 fim flux of Per-Bolo 58 is roughly consistent with that expected for the FHSC. The 24 /im detection is 
inconsistent with these simple spherical models of the FHSC, however. 
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Fig. 2. — Spitzer images of Per-Bolo 58 at 4.5,24,70 and 160 /an (left). Bolocam 1.1 mm continuum contours indicate the dense core, 
previously believed to be starless. The 70 /im map is from our new deep observations; other images are from c2d. At right is a 24 fim 
(grayscale) and 1.1 mm (contours) map of the surrou nding region in NGC 1333. 1.1 mm contours are 50, 100. ..200, 400. ..1400 mjy beam -1 
in both panels, and IRAS sources indicated are from Jennings et al. (1987). 
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Fig. 3. — Best-fitting FHSC model compared to the observed SED (left) and CARMA 2.9 mm visibilities (right; Schnee et al. 20TOD of 
Per-Bolo 58. While a PHSC provides a good fit to the SED, either a cavity in the inner envelope (R C ent ~ 50 AU), or outflow (a > 20deg) 
is required to match the 24 /im flux. The two best models, with Rcent = 50 AU, a = Odeg (red) and Rcent = 10 AU, a = 40deg (blue) 
are nearly indistinguishable. The input FHSC spectrum (L = 0.012 Lq, R = 1.7 AU, T e ff = 100 K; dashed line), and a model with no 
cavity (dotted line) are also shown for comparison. The millimeter visibilities help to constrain the envelope outer radius, with a best fit 
for R out = 8000 AU. 
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Fig. 4. — Same as Figure \3\ but for the best-fitting protostar model, with Lj nt = 0.012 Lq and a narrow outflow cavity (o ~ lOdeg). 
Despite the much hotter input spectrum (3000 K), the best-fitting protostar SEDs is nearly indistinguishable from the best-fitting FHSC, 
because most of the internal luminosity is re-processed by the dense envelope. Likewise, the 2.9 mm visibilities depend only on the envelope 
density profile, which is similar in both sets of models. 



